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FOREWORD 

This  report  was  prepared  for  the  Applied  Mathematics  Research  Laboratories 
by  P.  K.  Sen  and  P.  R.  Krishnaiah  under  Project  7071,  "Research  in  Applied 
Mathematics".  Part  of  the  work  of  P.  K.  Sen  was  performed  at  the  Aerospace 
Research  Laboratories  under  Contract  AF  33(615) -2272  of  the  above  laboratories 
with  the  University  of  Cincinnati.  The  work  of  P.  K.  Sen^in  the  final  stage, 
was  performed  under  Contract  F  33615-71-C-1927  of  the  Aerospace  Research 
Laboratories  with  the  University  of  North  Carolina. 
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ABSTRACT 


For  the  one-criterion  multivariate  analysis  of  covariance  (MANOCOVA)  model, 
the  rank  order  tests  for  the  overall  hypothesis  of  no  treatment  effect  considered 
by  Quade  (1967),  Puri  and  Sen(1969)  and  Sen  and  Puri(1970)  are  extended  here  to 
some  simultaneous  tests  for  various  component  hypotheses.  The  theory  is  based  on 
an  extension  of  rank  order  estimates  of  contrasts  in  multivariate  analysis  of 
variance  (MANOVA)  developed  by  Puri  and  Sen(1968)  to  the  MANOCOVA  problem,  and  is 
formulated  in  the  set  up  of  Gabriel  and  Sen  (1968)  and  Krishnaiah(1964 ,1969) . 
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1. INTRODUCTION 


Lee  Z^k)  =  (Y^k),X(k))'  =  (Y(k),...,Y(k),X(k),...,X(k))\  a  =  1 . n  be  n 

-a  -a  ~a  la  ’  pa  la  qa  ’  k  k 

independent  and  identically  distributed  random  vectors  (iidrv)  with  a  continuous 


cumulative  distribution  function  (cdf)  Il^(z)  ,  z  £  R^+C*,  the  (p+q) -dimensional 
Euclidean  space,  where  p  l,q  >_  1,  and  k=  l,...,c(  > 2 );  all  these  N  (= 


stochastic  vectors  are  assumed  to  be  independent.  The  q-variate  marginal  cdf  of 
X^k  is  denotea  by  G^(x),  x  £  R^,  and  the  p-variate  conditional  cdf  of  Y^k^ , 
given  X^k^  =  x  ,  is  denoted  by  F^(  y  Jx  ),  y  £  R*3 ,  x  £  R^ ,  k~l,...,c.  Our  basic 
model  is  the  following  : 


Gr(x)  =  G(x)  and  Fk(y  |  x)  -  F(  y  -  |  x)  ,  for  k=l,...,c,  (1.1) 

where  F  and  G  are  unknown  continuous  cdfs,  and  T,,...,  T  are  the  unknown 

~  1  -  c 

treatment  effects  (p-vectors) .  For  justification  of  this  model,  we  may  refer  to 
Scheffe^  (1959, ch. 6)  and  Sen  and  Puri(1970).  Tests  for  the  overall  hypothesis 


H  :  T  =  ...  =  T  =0  vs.  T,  j  0  for  at  least  one  k(  l<k<c)  ,  (1.2) 

o~l  ~  c  ~k  - 

based  on  suitable  rank  order  statistics,  were  considered  in  an  increasing  order 

of  generality  by  Qusde(1967),  Puri  and  Sen(1969),  and  Sen  and  Puri(1970).  We  are 

interested  here  in  testing  simultaneously  for  contrasts  among  t^,...,  t^,  both 

in  the  set  up  of  Gabriel  and  Sen(19X8)  and  Krishnaiah  (1964,1969).  These 

simultaneous  procedures  involve  the  concomitant  variates  adjusted  rank  order 

estimates  of  contrasts  among  t^,  which  are  considered  first  in  section  2. 

Section  3  deals  with  the  proposed  simultaneous  tests,  and  their  asymptotic  relative 

efficiency  (ARE)  results  are  presented  in  the  last  section. 


2.  COVARIATE  ADJUSTED  RANK  ORDER  ESTIMATES  OF  CONTRASTS  IN  MANOCOVA 
The  statistics  considered  by  Quade(1967),  Puri  and  Sen(1969)  and  Sen  and  Puri 
(1970)  ,  being  based  on  the  combined  sample  ranking,  are  not  suitable  for  our 
simultaneous  inference  procedures  to  be  considered  in  the  next  section.  For  this 
reason,  we  extend  the  results  of  Puri  and  Sen(1968)  on  rank  order  estimates  in 
MANOVA  to  the  MANOCOVA  problem,  which  provides  the  access  to  our  proposed  procedures. 
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Consider  the  pair  (k,£.)  of  samples,  and  let  n  ^  =  n^  +  n^  ,  for  1_<  k  <SL  <c. 

Let  R^k’^  be  the  rank  of  Y^k^  among  the  n.  0  observations  Y^k/ , . . .  .Y^  . 

i,a  ia  k£.  il  in^  ll 

,  for  a  =  1, . .  • ,  n,  ,  1  <  i  <  p ,  and  let  sfk’^  be  the  rank  of  X^k^  among  the 

in.  k  —  —  r  i,a  ia  6 

1  fW')  fk)  (1)  ( i' i 

n,  „  observations  X. . , . . . ,X.  ,  X..  , .  •  •  >  X.  -  a  =  1, . . . ,n,  ,  i=l . .  ,  for  1<  k 

kx  li  in,  il.  in.  k 

k  i 

<  l<c.  Consider  now  (p+o)  sets  of  rank  scores 

J^(a/(nk£+  1)),  a  =  1,...,  n  ^  ,  i=l, _ ,p  (primary  variate  scores)  (2.1) 

k£ 

*(i) 

(a/(n^^+  1))  a  =  1,...,  n^  ,  j=l,...,q  (  covariate  scores),  (2.^.) 

defined  for  every  (k,Sl):  1  <  k  <  SL  <_  c,  in  the  same  fashion  as  in  Puri  and  Sen 
(1971,  section  3.6).  Define  then  the  statistics 

Tu  ■  \  Cl  »>■  1 . p;  i;  K  <  1  i  c,  (2.3) 

kit 

TwJ)  ■  C  Cl  J^i)<sni*)/(v+l»-  i-1— (2.4) 

If  instead  of  yfk^  and  Y^  ,  one  works  with  yfk^  -  a,  a=l,...,n  and  Y^,  a  =  1, 
ia  ia  ia  k  ia 

...,n^,  the  corresponding  rank  order  statistic,  defined  by  (2.3)  is  denoted  by 

T^P(a),  1  <k  <  SL  <  c,  i=l,...,p,  and  <  a  <  00  .  We  now  assume  that  (u)  is 

nkX, 

t  in  u  :  0  <  u  <1,  so  that 


T^jp(a)  is  4-  in  a 


<  a  <  00 


(2.5) 


Nert,  we  rank  the  observations  coordinatewise  within  each  sample  separately.  Let 

R*(k)  (  or  S*<k>  )  be  the  rank  of  Y^k)  among  Y^k)  , . . . ,yfk)  (  or  xfk)  amnong  X(k) , . 
i,a  i,a  ia  6  il  in.  ia  il 

(k)  K 

X)  ),  for  a  *  l,...,n  ,  i*  (or  i-  and  k*  Define 

|C 

,  1*1,..., p,  and  j=l,...,q  as  in  (2.1)  and  (2.2)  with  the  n  .  being 

nk  "k 

replaced  by  n^  ,  and  let 


v*k)  =  nf1  J(i)(R*(k)(n.+l))J(:1)(R*(k)(n,+l))  -  J(i)J(j),  (2. 


'ij  "k  a=l  Jnk  i,a  k  \  '  j.cf'  k 
for  i,  j  =  1,...,  p,  where 


nk  nk 


6) 


T(i)  _  r  k  t(D 


n.  I  *  Jv  '(a/(n  +  1)),  i  =  l,...,p. 

1C  Ct  i  n,  k 

k 


(2.7) 


Similarly,  define  for  p+i  ^  i,  j  p+q  by  replacing  J^'*  and  R.^  by  j 

*CW>  ^  ^  ’  k 

and  S.  '  respectively  (  and  also  for  j)  in  the  definitions  in  (2.6)  and  (2.7). 
1 ,  ot 

Finally,  keeping  and  R.  ^  as  they  are  while  replacing  and  R.^  by 

nu  1  •  ni,  J  »a 

J  ''■''and  S.  .respectively,  in  (2.6)  and  (2.7),  we  define  in  the  same  way  v.. 
nk  j,a  ij 

for  i=l,...,p;  j=p+l, . . . ,p+q.  Let  then 

V(k)  =  ((  vjk))).  .  ,  k=l,...,c,  (2.8) 

~  lj  1  ,  J_l,  •  •  •  .P"™! 

( k )  fk') 

where  v.  ,  =  v.  /  for  all  i,j.  Further,  let 

Ji  ij 


h  -  £k.iVN)r 


/  ~N,11  ~N,12  \ 
V  ~N,21  ~N,22  / 


where  V  , ,  ,  V. ,  .  „  =  V  and  V„  are  respectively  of  the  order  pxp,  pxq  and 

^<N )  il  ^>N  j  1Z  ,  l. x 

qxq.  Finally,  let 


V  =  V  -  V  V  V 
IN  -N,  11  -N,12~N ,22  4N,21  ’ 


(2.10) 


where  V  ??  is  a  generalized  inverse  of  V  99  .  Define  then  J^1  ,  i=  l,...,p,  and 

— *(  i)  KX# 

J  J  ,  j=l,..., q,  for  1  <_  k  <  £  c  as  in  (2 . 7)  with  n  being  replaced  by  n  ,  , 

nk«,  k  kx- 

and  let 


‘  <  . Tu<V  >•  •  *5  k  <  11  i  c  ; 


-  <  T. 


, ,  x*^  )'j  1  <  k  <  £  <  c; 


(2.11) 


(2.12) 


»  <  jJP . 'if’)’.  3?  -  (  li11 . ffl'.l'l'l  £  c;  (2.13) 


&<•>  ■  <  Ik*<?>  -  »  -  X»,12Y»;22(  Ik*  -  tl  > 

'  <  Tk*1)(*l>  . TUP>(V  >’•  f"  1ik  <  1  -  c. 


(2.14) 


The  statistics  in  (2.14)  are  the  covariate  adjusted  rank  order  statistics  as  defined 
in  Sen  and  Puri(1970)  with  the  only  difference  that  their  combined  sample  permuta¬ 
tion  covariance  matrix  is  replaced  here  by  V^,  the  average  within  sample  rank  order 
covariance  matrix.  This  replacement  is  made  with  the  primary  objective  of  making 


the  adjustments  for  the  covariates  in  (2.14)  unaffected  by  t,,...,  is 


translation  invariant,  so  that  the  unknown  T ^ . have  no  effect  on  it,  while 

the  combined  sample  permutation  covariance  matrix  rests  on  the  assumption  that 

1 


x,  =  ...  =  t  ,  which  is  not  necessarily  true. 


~c 


Now,  precisely  by  the  same  alignment  logic  as  in  Puri  and  5en(1968) ,  we  define 

*(i) 


k£ 


=  inf {  a:  T°(l)(a)  <  0  },  A^^  =  sup  {a:  T°(l)(a)  >  0  },  (2. 


15) 


for  1  _<  k  <  £  <c,  i=l,...,p.  Then,  our  proposed  estimator  of  ^ 


=  T, 


iit  ■  <"4k“ . *w' >'=  Ki  ■  t  <  ‘i 1 ;  »• 


(2.16) 


Conventionally,  we  let  A 


kk 


Akk  =  k=1»---»c  and 


-  tor  1  <k  <  L  <c. 


Then,  as  in  Puri  and  Sen  (1968;  1971, Ch.  6),  we  define  the  compatible  estimators 
'* 


Ak£  =  -k.  -£.  5  -k.  c  Z£=l  -k£  ’  1  -  k  *  1  -  c’ 


(2.17) 


Let  H,..(x)  (  and  G,..(x)  )  be  the  marginal  cdf  of  Y^  _  (  and  xfk^  ), 

(l)  (j)  tot  k  ia 

i=l,...,p  (  and  j=l,...,q),  and  let  H^^.^x.y),  H^^(x,y)  and  G(j  ^(x.y) 

be  respectively  the  joint  cdf  of  (Y^k)-  T^l)  ,  Y^k)-  '  ^ )  ,  (YW-  T^l)  ,  X^k)  ) 

r  ia  k  i  a  k  ia  k  ja 

and  (X^k\  X<«  ),  for  i(  ¥  i')  =  l,...,p  and  j(  ^  j')  =  l,...,q.  Define  then 

J(l)(u)  =  Jll)(u>  ’  J*(j)(u)  =  Jn(3)(u)  ’  0  <  u  <  M2. 18) 

n-xjo  n  n^°°  n 


for  i=l,...,p,  j=l,...,q,  and  denote  by 

*  (i)  *  *(i) 

=  /„  JV  '(u)du  ,  i=l,..,p;  Uj  =  / c  J  vj;(u)du  ,  j=l,...,q; 

“ii.u  ■  £  £  Jt«(H(1)(,))J«J!<»(J|W)d.«)(,,,) 

-  vt  Wj  ,  M  =  ; 

,12  =  C  C  J(i)(H(.)(x))J*(j)(G(.)(y))dH^>.)(x,y) 

* 

-  uiPj  ,  i=l,...,p;  j=l,...,q; 

Vij,22  =  C  C  J*(i)(G(i)(x))J*(j)(G{j)(y))dG(i>j)(x,y) 


V.  . 

1J 


('u 

*12 

V  = 

'  -21 

-22 

) 


-  u*i  dj*  •  M=  M..,q; 


-11  "  ((vij,ll))’  -22  ((vij,22)) 


(2.19) 


(2.20) 


(2.21) 

(2.22) 

(2.23) 
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and  v  =  v  =  ((v..  )).  Further,  let 

~12  ~21  ij  > 12 

*  —  ★ 

V  =  V  -  v  v  =  ((v..)).  .  ; 

~  ~11  ~12  ~22  ~21  ij  i,j=l,...,p 

r*  =  ((  Y* .) )  =  ((  v.  ,/B.B.)).  . 

ij  ij  i  1  i.J=l.- • •  .P 

where 

,oo  ( t  } 

B.  =  f-ao  (d/dx) J'  ;(H(i)(x))dH(i)(x)  ,  i  =  1.....P. 
Finally,  consider  a  contrast 

<f>(£)  =  2.x.  +  ...  +  £  t  ,  £  =  (£ . %  1  1  =  (1,. 

~~  i~i  c~c~  1  c  ~ 


.,1)’, 


(2.24) 

(2.25) 

(2.26) 

(2.27) 


for  which  we  have  the  compatible  estimator 

* 


<t>  (£)  =  £,  A,  +...+£  A  ,  (2.28) 

where  the  A^  are  defined  by  (2.17).  For  any  £  J_  1,  (2.28)  provides  a  robust, 

traslation  invariant  and  concomitant  variate  adjusted  estimator. 

* 

For  the  study  of  the  asymptotic  properties  of  <J>  (£)  ,  we  assume  that  as  N  ■*  °°, 
n^/N  -*■  :  0  <  <  1,  for  k  =  l,...,c.  (2.29) 

Then,  by  the  same  technique  as  in  theorems  2.1  and  2.2  of  Puri  and  Sen(1968)  (  with 
direct  adaptations  from  theorem  5.1  of  Puri  and  Sen(1966) ,  yielding  the  joint 
asymptotic  normality  of  the  covariate  adjusted  rank  order  statistics  in  (2.14)), 
we  obtain  that 


X#(  N1/2[  $*<£)  -  $(£)]  )  -  X(  0,  (  (2-30) 


Similarly,  if  we  have  r(  >  1)  linearly  independent  contrasts 


-  \=i  *sk  ik> 

£ 

-s 

1  1  ,  s  =  1, — ,r  ,  (2.31) 

and  we  define  C  =  ((  c  ,))  , 

SS  S  y  S 

. .  ,r  by  Css'  =  Ek=l  £sk£s’kAk’  s-s'=1»-*-»r. 

we  have 

£,(N1/2[  Ay  -  My],  1  <s<  r)  -*■  ^(D,  C8  r*),  (2.32) 


where  fi  stands  for  the  Kronecker  product  of  two  matrices. 

We  shall  make  use  of  (2.30)  and  (2.32)  for  providing  suitable  confidence  inter- 
vals  to  <£(£)  ,  £  X  1  .  For  this,  we  require  to  estimate  the  unknown  F  .  Now,  on 


6 


defining  t  as  the  upper  100a%  point  of  the  standard  normal  distribution,  and  letting 


~  sup*  a:  TU  (a)  >  tn/2[  VN,ii(nk«,/nknP ]  }  * 

Ak£,L  =  lnf{  a:  Tk  (a)  <  ta/2f  VN,ii(nk«/nkn£^  }  ’ 


(2.33) 

(2.34) 


and  following  the  same  method  of  proof  as  in  theorem  1  of  Sen(1966),  it  follows  that 

■  (2ta/2<\.u)1/2)/[(V/V1l)1/2<  ~AU,U  ~  <2'35> 

are  all  translation  invariant  consistent  estimators  of  B^,  defined  by  (2.26),  for 
i=l,...,p.  Hence, 


(S-1  Z  B(i) 

'2'  Kk<l<  c  k£ 


,  i  1 , .  •  •  ,  P 


(2.36) 


are  translation  invariant  robust  and  consistent  estimators  of  B^,...,B^,  iespectively . 

Also,  by  theorem  4.2  of  Puri  and  3en(1966)  [  see  also  theorem  3.2  of  Ghosh  and 

Sen(197])]  and  some  standard  manipulations,  it  follows  that  whatever  be  I, . T  , 

~  1  ~  c 


v  as  N  ■*  oo  , 


(2.37) 


Thus,  from  (2.10),  (2.24),  (2.25),  (2.36)  and  (2.37),  we  obtain  on  defining 


rN  =  ((  Yt1  »  ;  Ytl  =  V  /[  B(1)B(J)],  i.j  =  1 . . 


(2.38) 


as  N  +  ” 


(2.39) 


Now,  the  space  of  all  possible  contrasts  $(£)>  £  J_  1,  is  spanned  by  a  set  of 
c-1  linearly  independent  contrasts.  Hence,  using  (2.31)  and  (2.32)  with  r  =  c-1,  and 
(2.39),  we  obtain  by  Roy-Bose  method(1953)  that 


lim  p{  Nl/2f  t,(  _  <{>(£))][ (t ' ?*t)  (Ek=1^/^k>  ] ”1/2 


(2.40) 


—  ^p(c-l),a  ^°r  evety  ~  J  ^  9  ^  1-a, 


which  provide  a  simultaneous  confidence  region  to  the  set  of  all  possible  contrasts 
2 

where  \  is  the  upper  100  %  point  of  the  chi-square  distribution  with  p(c-l) 

P  \  C  L  )  9  Ct 

degrees  of  freedom.  ARE  results  will  be  considered  in  section  4. 
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3.  SIMULTANEOUS  TEST  PROCEDURES 

Along  the  lines  of  Krishnaiah  (1964,  1969),  we  consider  here  some  simultaneous 
test  procedures  (STP)  for  contrasts  in  MANOCOVA  based  on  the  robust  estimators 
derived  in  the  earlier  section. 

Procedure  I.  Consider  a  set  of  r(>  1)  linearly  independent  contrasts  (vectors) 
$(^s) ,  s=  l,...,r,  defined  as  in  (2.31),  and  let 


♦  u8>  =  0  vs  As:  4><y  ^  0  >  s  =  l,...,r  ; 


(3.1) 


H  =  H1  Hr  and  A  =  Ax  O  . . .  (J  Ar  . 


(3.2) 


Let  us  define  the  compatible  estimators  <j>  (f.g) ,  s  =  l,...,r,  as  in  (2.28),  and  let 


Vs  ■  • s- 1 . r- 


(3.3) 


where  is  defined  by  (2.38),  and  csg»  in  between  (2.31)  and  (2.32).  By  using 
(2.32),  (2.39)  and  a  result  of  Krishnomoorthy  and  Parthasarathy  C1951),  we 
conclude  that  under  H  in  (3.2),  Qjj  =  ^  Qjj  !»•••»%  r)  '  has  asymptotically  a 


multivariate  chi-square  distribution.  Hence,  there  exists  a  such  that 


lim 

N-WO 


p{  Qn,s  1  Qa  >  Vs  =  I  H}  =  1  -  a 


(3.4) 


Our  proposed  STP  consists  in  rejecting  Hg  in  favor  of  Ag  only  for  those  s  for 

which  QN  s  >  Qa  ;  otherwise,  accept  Hg,  s  =  l,...,r.  The  total  hypothesis  H  is 

accepted  iff  all  the  component  hypotheses  are  accepted.  The  associated  simultaneous 

confidence  intervals  for  t'<J>(£  ),  t  f  0,  are 

-  ~  ~s 

|t'(  ♦*«  )  -  «K4B))  I  <  {Q  c  (t'r*t)/N  }1/2,  s  =  1, . . .  ,r,  t  i  0,  (3.5) 

with  an  overall  confidence  coefficient,  asymptotically ,  equal  to  1  -  a  . 

Procedure  II.  This  is  really  a  stepdown  procedure.  We  denote  by  CL  the  principal 

minor  of  C  comprising  of  the  first  j  rows  and  columns,  and  write  for  2  £  j  <  r, 

fC.  ,  c. 


-i 


-  c’  ,C.  ,c.  ,,  j  =  2,...,r,  and  let  c 


11  =  C11  ‘ 


Then,  by  (2.32),  on  denoting  by  =  N^2[$  (£^)  -  4>(£^)]  ,  and  for  j  ^  2, 
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,1/2 


z  =  (a1>  -  <$  (Vj.j)  -  ^aj_1)),],> 


(3.6) 


we  conclude  that  has  asymptotically  a  multinormal  distribution  with  mean  0 

* 

and  dispersion  matrix  c, ,  while  given  Z.,  the  conditional  distribution  of 

11  ~N  -3 

1/2  * 

N  [  4>  (Z  )  -  ]  is  asymptotically  multinonral  with  mean  vector 


and  dispersion  matrix 


*  'k 

c.  .  f.,  ,  for  j  =  2,. . .  ,r. 


jj  ~N 


(3.7) 


(3.8) 


Let  us  then  write  =  $(£^)  and  for  2  <  j  r, 

n.  =  $(Z  )  -  (cj_i  C"^  ®I)[($(Z1)),,.-.,(*(Z  x))']’. 

Then,  the  total  hypothesis  H  in  (3.2)  may  be  written  as  H  =  (1 , 

* 

:  Hj  =  0  ,  for  j  =  l,...,r. 


(3.9) 


O  H  ,  where 
r 


(3.10) 


We  cfesire  to  provide  a  STP  for  H.  ,  j=  i,...,r.  For  this,  let 

*  * 


»  I1 


hj  =  4>  (Z  )  -  (cj^  C  &I  )[($(Z1))',...,  (♦(Zj_1))'] 


(3.11) 


for  j  =  2,....,r  and  = 


(Z^) ,  and  let 


QN,j  =  (N/C3*3H(  ~j),(  1(  Dj)]  •  J  =  1 . r- 


(3.12) 


By  (2.32),  (2.39),  (3.7)  and  (3.8),  we  conclude  that  under  H,  Q, 


N,1 


QN>r  are 


asymptotically  independent,  each  having  a  central  chi-square  distribution  with 


p  degrees  of  freedom.  Hence,  there  exists  a  such  that 

p  {  V.  i  <  ■ !  ■ 1 . 'I  H 1  • 1  -  “  • 


(3.13) 


*  * 

Our  proposed  STP  consists  in  testing  for  , . . . ,  sequentially.  That  is,  if 

*  *  *  *  * 

Q„  ,  >  Q  ,  reject  H,  and  hence  H  ;  otherwise  proceed  to  test  for  H„.  If  0„  „  > 
N,1  a  J  1  2  N,2 

*  *  * 

,  reject  and  hence  H;  otherwise  proceed  to  ,  and  so  on.  The  total  hypo- 

*  * 

thesis  H  rs  accepted  iff  for  all  j  =  l,...,r.  By  virtue  of  (3.7)  and 

(3.8),  we  obtain  the  following  associated  simultaneous  confidence  intervals 


|n1//2(  n*  -  o .)  't 

~J  ~3  ~ 


<  [  Q*  (  t'(?*)t  )c..]1/2  j  =  1 . .  t  t  0, (3.14) 


•a  '  ~  VN  ~ 

with  an  asymptotic  confidence  coefficient  equal  to  1  -  a 
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Procedure  III.  Let  x'  =  (  x'  .  x')  .  and  let 

-  -  ~  I  ~c 


(  an,  a12,...,alp,  ....  acl,  acp)  +  O' 


(3.15) 


be  a  pc-vector  such  that  a^,...,  a^  are  not  all  equal  for  1  =  l,...,p;  the 
totality  of  all  possible  a  satisfying  the  above  condition  is  denoted  by  A  .  Then, 
the  hypothesis  H  in  (3.2)  may  be  written  as  H  =  H  ,  where 


Ha  :  a'x  =  0  for  a  e 


We  want  to  consider  a  STP  for  all  H  ,  a 

a  ~ 


(3.16) 


tA»  • 


Let  us  denote  by  A*  =  (  A'  ,  ....  A'  ),  where  the  A  are  defined  by  (2.17), 

~  .L  •  ~  C  •  ~  K.  • 

and  let  n  =  Diag(  n^,...,  nc) •  Define  then 

Vo  ■  ®<Srl>  l  -  k.&'L-  °-in 

By  (2.32)  and  (2.39),  under  H,  Q„  has  asymptotically  chi-square  distribution  with 

N,o 

p(c-l)  degrees  of  freedom.  Hence, 


lim  P  {  Q 
N  -*■  00  TJ,o 


*-p(c-l)  ,a 


H  }  =  1  -  a 


(3.18) 


Replacing  A  by  A  -  x  in  (3.17),  and  thereby  eliminating  H  in  (3.18),  we  obtain 
by  the  Schwartz  inequality  that  asymptotically  with  a  confidence  coefficient  1  -  u, 

|a'(  A  -  x  )|  £  Xp(c_1)>0(  U(n  1(g)  T*  )  a  ]1/2,  \/  a  e  eA  •  (3.19) 

which  provides  a  simultaneous  confidence  region  to  all  possible  a'x  ,  a  e  • 

The  STP  consists  in  rejecting  those  H&  (  a  eA)  for  which  | a* A |  >  Xp(c_i)>a* 

[a'(n  ^^r^)a]^2.  A  similar  procedure  for  MANOVA  (  but  involving  a  loss  of  A.R.E. 
for  contrasts  not  involving  all  the  c  samples)  is  due  to  Gabriel  and  Sen(1968) . 

Procedure  III  has  the  maximum  flexibility  to  all  a  6  JL  ,  but,  in  practice, 
when  one  may  be  interested  in  a  set  of  specified  contrasts,  it  is  usually  less 

* 

efficient  than  the  two  other  procedures.  The  simplicity  of  the  computation  of 

*  * 

over  Q  (  derived  from  the  asymptotic  independence  of  Q  QN  )  has  to  be 

weighed  against  the  arbitrariness  and  complications  invo  .ved  in  the  choice  of  the 

£  A 

sequence  of  the  suffixes  l,...,r  and  the  formulation  of  One  may  also 

consider  some  other  procedures  along  the  line,?  of  Krishnaiah(1969) . 
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4.  ARE  RESULTS 

One  could  have  Ignored  the  covariates  totally,  and  proceeding  as  in  Puri  and 
Sen(1968) ,  obtained  rank  order  estimates  of  the  contrasts  derived  solely  from  the 

1  ^  k  <  £  £  c,  i  =  l,...,p.  In  this  case,  the  asymptotic  covariance  matrix 

* 

in  (2.30)  or  (2.32)  would  have  involved  T  instead  of  T  ,  where 

r  =  ((  Y  ))  ,  i.j  =  !»•••. P.  (4.1) 

and  V„  n  a  ^  are  defined  by  (2.20)  and  (2.26).  Thus,  as  in  Puri  and  Sen(1968) , 
defining  the  .ARE  as  the  reciprocal  of  the  ratio  of  the  generalized  variances  raised 
to  the  power  p  \  the  ARE  of  the  covariate  adjusted  rank  order  estimates  with 
respect  to  the  unadjusted  ones  is 


«!  -  <l[|/|r*|}1/p  -  <|vu|/|v*|  }' 


^-ll^-ll  "  -12-22-21  P’  (4*2) 

where  V^2V22J21  is  Pos:'-tive  semi-definite.  Hence,  we  have 

e^  _>  1,  where  the  equality  sign  holds  when  =  0  •  (4.3) 

This  explains  the  asymptotic  supremacy  of  the  covariate  adjusted  estimates.  It  also 
follows  similarly  that  the  ARE  is  a  non-decreasing  function  of  the  number  of 
concomitant  variates,  i.e.,  additional  covariates  induce  more  information  on  the 
estimates. 


For  the  classical  normal  theory  estimates  of  contrasts  based  on  the  sample 


mean  vectors,  let  Z 


and  let 


hi  -12 
hi  -22 


be  the  covariance  matrix  of  (Y^  '  ')  '  , 

-a  -a 


l  -11  hlhllll  '  (4,4) 

Then,  from  the  results  of  Sen  and  Puri (1970),  we  may  conclude  (  on  omitting  the 

* 

details)  that  the  ARE  of  <J>  (£)  with  respect  to  the  classical  normal  theory  estimator 


of  <K&)  is 


{  |i*|/|r*|}1/p, 


which  depends  on  ^n’  ^22’  ^12  aS  we^  as  ^ll’  -22’  V12  an<^  B^,...,B  .  In 
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general,  it  is  not  possible  to  attach  suitable  bounds  to  However, 

{|E*l/|r|)1/p  <  e2  <  {|En|/|r*|)1/p.  (4.6) 

In  particular, if  £  =  0,  e„  >_  {|L|  |r|}^p,  where  the  equality  sign  holds  if 

v^2  =  0  •  Thus,  if  the  original  variates  are  uncorrelated  but  not  necessarily 

independent,  we  expect  to  gain  in  efficiency  by  considering  the  robust  estimates 
* 

<j>  (8-).  The  same  ARE  results  holds  for  the  simultaneous  confidence  intervals  in 
(2.40)  when  compared  with  the  parallel  procedure  for  the  MANOVA  model  and  the 
parametric  Roy -Bose  type  procedure. 

For  the  simultaneous  tests  in  sectiou  3,  the  ARE  with  respect  to  their 
parametric  counterparts  considered  in  detail  in  Krishnaiah(1969)  depends  on  the 
particular  set  of  variates  or  samples  included  in  the  set  of  contrasts  under  test. 

The  overall  ARE  agrees  with  e^  or  e^  in  the  respective  cases,  while  the  minimum 

and  maximum  (  over  all  possible  choice  of  variates  and  samples  )  ARE  are  respectively 

*  * 

the  minimum  and  the  maximum  characteristic  roots  of  E  (  T  )  .  Since,  these 

bounds  are  similar  to  the  ones  discussed  in  Sen  and  Puri  (1970) ,  we  omit  the  details 
here.  In  passing,  we  may  remark  that  on  using  the  normal  scores  statistics  for 
deriving  the  robust  estimates  in  section  2,  we  are  able  to  achieve  asymptotic 
optimality  when  the  underlying  distribution  is  normal. 
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